Limits to phase resolution in matter wave interferometry. 
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We study the quantum dynamics of a two-mode Bose-Einstein condensate in a time-dependent 
symmetric double-well potential using analytical and numerical methods. The effects of internal 
degrees of freedom on the visibility of interference fringes during a stage of ballistic expansion are 
investigated varying particle number, nonlinear interaction sign and strength as well as tunneling 
coupling. Expressions for the phase resolution are derived and the possible enhancement due to 
squeezing is discussed. In particular, the role of the superfluid - Mott insulator cross-over and its 
analog for attractive interactions is recognized. 



PACS numbers: 03.75.Be, 03.65.Ge, 05.60.Gg 

I. INTRODUCTION 

Coherent atom optics offers considerable promise for 
applications in a numbers of areas from precision mea- 
surements to rotation sensors, accelerometers, and grav- 
ity gradiometersP, |2- One key element in any practical 
device is a coherent beam splitter, and much effort has 
been devoted to the realization and the understanding of 
these devices, both experimentally and theoreti- 

cally illiiliaillllllll. 

In contrast to light fields, which do not interact in a 
vacuum, matter waves are subject to collisions, mostly 
two-body interactions in the low density beams normally 
considered when using quantum-degenerate atomic sys- 
tems. At extremely low temperatures, collisions produce 
a nonlinear phase shift of the matter waves that is pro- 
portional to the atomic density, and hence leads under 
normal conditions to undesirable phase noise This 
is a serious difficulty that needs to be addressed in detail. 
On the other hand, two-body collisions are also known 
to act as the matter- wave analog of a cubic nonlinearity 
in optics. As such, they can be used to generate nonclas- 
sical states of the Schrodinger field. These states can in 
turn be exploited to achieve phase resolution below the 
standard shot-noise limit. 

The use of squeezing to reduce quantum noise was sug- 
gested in [Tsj l for Ramsey-type interferometry. For optical 
Mach-Zender interferometry a Heisenberg limited scheme 
was outlined in (16] , and related schemes using dual Fock- 
states were suggested for Bose-condensed atoms later. 
fvA ITM IT^ . A scheme dependent on parity measurement 
was suggested in poj . 

A matter-wave beam splitter can be thought of as a 
double-well potential with time-dependent well separa- 
tion. At zero temperature, the dominant mechanism by 
which atoms move from one well to the other is quan- 
tum tunnelling. It is known that in quantum-degenerate 
bosonic systems, the interplay between tunnelling and 
collisions and the associated mean-field energy can re- 
sult in highly non-trivial effects 21] . For instance, in the 
case of repulsive interactions a condensate trapped on a 
lattice potential can undergo a quantum phase transi- 



tion from a superfluid state, characteristic of situations 
where tunneling is dominant, to a Mott insulator state, 
characteristic of situations where the mean-field energy 
dominates the dynamics I23L |23 ]. The situation is 
more complicated in the case of attractive interactions, 
since the condensate is then unstable unless it is small 
enough to be stabilized by its kinetic energy |2^. Such 
small, stable condensates also undergo a transition remi- 
niscent of the superfluid-Mott insulator transition in the 
sense that the ground state ^EE3 changes its statistical 
properties at a critical value of the interaction strength. 

With these processes in mind, the goal of this paper 
is to assess in detail the limits in phase resolution of 
an atomic beam splitter under the combined effects of 
collision-driven cubic nonlinearities and quantum tun- 
neling. The combination of these effects renders a full 
quantum-mechanical description of an atomic beam split- 
ter highly non-trivial. Consequently, we restrict our dis- 
cussion to a two-mode analysis, using a combination of 
numerical and analytical tools. 

This paper is organized as follows: Section UTI discusses 
our model and establishes the notation. Section IIII Al 
presents results of a numerical analysis of the static prob- 
lem where the condensates are released from a beam split- 
ter with fixed well separation. Depending upon the ratio 
of mean field energy to inter-well tunneling energy, the 
beam splitter operates either in the superfluid or Mott 
insulator-like regime, with qualitative and quantitative 
differences in their noise properties. For attractive in- 
teractions, the phase noise of the beam splitter is found 
to be significantly reduced at the transition between the 
two regimes. These results are extended to the dynami- 
cal regime in Section IlII Bl which discusses in particular 
the departure of the system from adiabaticity. Finally, 
section Hvl is a summary and outlook. 



II. MODEL 

We consider in H-l dimensions the quantum dynam- 
ics of an ultracold bosonic atomic beam trapped in a 
double-well potential V{y, d) with time-dependent well 
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separation 2d{t). A beam splitter using such a time- 
dependent configuration of optical waveguides has been 
realized by the MIT group [23 , the resulting atomic field 
being detected after turning off the trap and ballistic free 
expansion of the atomic condensate. Assuming that the 
atomic density is low enough that wc can neglect three- 
body collisions, the Hamiltonian of this system is 
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dy¥{y)¥{y)i>{y)i>{y), 
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where 52 is the two-body coupling constant, taken to be 
negative for attractive two-body interactions. In the ab- 
sence of excitations to higher spatial modes, the field 
operator can be expanded into two modes corre- 

sponding to particles located around the two minima of 
the double-well potential as 



= VL{y, d)dL + (pniy, d)aR. 



(2) 



Here di(fl) are bosonic annihilation operators for the 
"left" and "right" mode of the matter-wave field, and 
fL{R) are the corresponding spatial mode functions. For 
sufficiently harmonic potentials they can be approxi- 
mated by 



fiv, d) 



1 



L{R) 



(irAy^) 4 



cxp 



4Ay2 



(3) 



the minus sign corresponding to ipii{y,d) and the plus 
sign to (fLiy,d). Within this two- mode approximation, 
the Hamiltonian becomes 



H{t) = huj{d\dL + a^i^aii) + ^-^^^^ [a[a_B, + o^Ql] 
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where we have introduced the time-dependent tunnelling 
energy 



AEit)^ j dy^L{y,d{t)) + V {y,dit))^ ipRiy,dit)) = hu;exp[-d\t)/Ay% (5) 
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and g — (72*^4 with 

d4- / dyipUy,dit)) = -^—. (6) 

Note that we have neglected cross-phase modulation, 
consistently with the validity of the gaussian approxi- 
mation in the description of the mode functions of the 
waveguide. We remark that this approximation only 
holds for d > Ax, otherwise the modes must be taken as 
time-dependent linear combinations of the energy eigen- 
states, a procedure requiring numerical diagonalization 

0- . 

This two-mode problem is conveniently reexpressed 
in the Schwinger angular momentum representation of 
bosonic operators 29]. We proceed by introducing the 
angular momentum operators 

Jz = ^{aldL~d\^dR), (7) 

= ^{dlcLR + cl^cil), (9) 

which can be thought of as the orthogonal components 
of a Bloch vector of length N/2. This corresponds to 



mapping the quantum state onto a distribution on the 
Bloch sphere. 

As usual, we then express the state of the matter- wave 
field in terms of eigenstates | J, m) of the operators ,P 
and Jz, where 

i2 = Jl + Jl + Jl (10) 

with 

P|J,m) = h^J{J+l)\J,m) 

Jz\J,m) = hm\J,m) ^ ^{hl - nR)\J,m), (11) 

and J = N/2, m = - J, - J + 1, . . . , J. 

In the angular momentum representation the Hamilto- 
nian Q reads 

f{J) + 2gJl + AE{t)J,, (12) 

where the energy /(J) is a function of the total angular 
momentum eigenvalue J. For a fixed particle number, it 
yields a constant phase shift irrelevant for the problem 
at hand. The ground state of the Hamiltonian (|12|l is 
expressed in terms of the azimuthal quantum number 

J 

Cm|>/,m). (13) 

m— — J 
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Each of the operators 0-® generates rotations of this 
distribution around the corresponding axis. As seen from 
the Hamiltonian (|12|l . and aheady proposed in Ref. 
a rotation about the x axis of the Bloch sphere can be 
achieved by turning on the quantum tunnelhng between 
the two wells for a precisely determined time. As a re- 
sult, it is possible to transform a number-squeezed state, 
characterized by reduced fluctuations in ,into a phase- 
sc(ueezed state, characterized by reduced fluctuations in 
Jy. We exploit this feature of quantum tunnelling later 
on in to achieve sub-shot noise detection in the presence 
of repulsive interactions. 

We mentioned that the detection of the atomic field is 
carried out after the optical waveguide is rapidly switched 
off and the atoms undergo a stage of ballistic expansion. 
The mode functions ^ no longer describe the spatial 
density of the condensate during that stage. Rather, 



they must be replaced by free gaussians that are cen- 
tered around the minima of the potential at the time of 
release 

X -P (- ^t'^' . (14) 

Taking the two halves of the condensate to have a relative 
phase G, the field operator for the ballistically expanding 
atoms becomes then 

i!{y,t,Q) = aL¥>L(y,t)exp(ie/2) 

+ ai?¥>fl(2/,Oexp(-ie/2), (15) 

resulting in the spatial density 
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The atomic density at any point in space and time is 
given by an incoherent contribution that is independent 
of both the relative phase and the internal dynamics of 
the two-mode condensate, as well as a coherent contribu- 
tion. 

Since the Hamiltonian H is invariant with respect to 
the exchange L ^ R, it is easily seen that for states sym- 
metric with respect to the interchange of the two modes, 
(Jy) = (J^) = for all times, although (J^) ^ as we dis- 
cuss shortly. The coherent contribution to {Gi{y, t, 8)) is 
therefore proportional to the expectation value of (Jx)- 
In terms of the angular momentum picture, it can be 
interpreted as the polarization of the distribution. The 
visibility of the interference fringes of the ballistically ex- 
panding two-mode condensate, 



V = 



(Gi(0,^,0))-(Gi(0,t,7r)) _ 1(4)1 
(Gi(0,i,0)) + (Gi(0,i,^)) A^/2 ' 
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depends only on that expectation value, the associated 
fluctuations being given by 



AGi ^ ^{AGl{y,t,Q)) 



= ^{Gl{y,t,e))^{G,iy,t,Q)y- (18) 

The phase resolution of atom-interferometric experi- 
ments is limited by the requirement that the change in 
local density resulting from an imprinted global phase 
change must be larger than the intrinsic fluctuations of 



the first order correlation function 

d{G, 



AGi > 



AO. 



(19) 



This gives as an estimate for the phase resolution result- 
ing from number fluctuations 



Ae = 



V(AGf) 

|9(Gi)/ae|- 



(20) 



Applying this criterion to the problem at hand, one 
finds that the phase resolution for the two-mode conden- 
sate is given by 



Ae^ 



(j2) , {N^)-{Ny^ 



(21) 



The second term in Eq. (|21l) is equal to zero in an indi- 
vidual experimental event, but gives rise to a contribution 
of order unity if averaged over a classical (Poissonian) dis- 
tribution of particle numbers in repeated experiments. 
The number of particles must therefore be determined 
with a high precision if one wishes to benefit from any 
enhancement of interferometric resolution due to phase 
squeezing. Setting this second term to zero for now, we 
recover the result derived by Kitagawa and Ueda in their 
seminal work on spin squeezing |3ll | . Squeezing along the 
y-direction corresponds to increased correlations of the 
phases between the two wells, thus producing a better 



4 




FIG. 1; Particle density bounded by density noise (Gi)±AG'i 
after ballistic expansion of a two-mode condensate with TV = 
100 in the 'superfluid' regime where the visibility is high. The 
dashed lines show the density of the mode functions. Trans- 
verse distance is measured in units of the oscillator width Ay 
and the initial condensates were centered around d = ±8Aj/. 



defined relative phase, a desirable feature in interfero- 
metric applications. The squeezing parameter, which we 
calculate relative to the standard quantum limit 
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although not necessarily representing the maximal 
achievable squeezing of the variances of the angular mo- 
mentum distribution, was shown by these authors to rep- 
resent the parameter of interest from the point of view 
of possible interferometric applications. Hence it is the 
major focus of the present study. 



III. RESULTS 

This section discusses the main results of our numerical 
study of the two-mode beam splitter. The analysis is 
based on the dimcnsionless parameter 
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2gN 
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the ratio of mean-field energy to tunnelling energy, which 
completely determines the spectrum of the Hamiltonian. 

Note however that the two-mode approximation im- 
plicitly assumes a stable condensate, a property that 
holds only for low enough densities in the case of at- 
tractive interactions p5l |. 

We consider first the simple case of a static double- 
well system where the pconfining potential is suddenly 
turned-off, resulting into the formation of an interference 



pattern after a period of ballistic expansion. The next 
subsection discusses the results of a full dynamical study. 



A. Static double-well potential 

The matter-wave interference pattern is shown in Fig. 
^ in the superfluid regime, further illustrating the ex- 
cellent contrast in that case. In addition to the fringe 
contrast, it is necessary to consider the quantum fluctu- 
ations of the interference pattern, since they lead to the 
fundamental limit in phase resolution of interferometric 
measurements. Figure ^ illustrates these fluctuations by 
attaching to the intensity (Gi) a width given by twice its 
variance AGi = ((G^) — (Gi)^)^/^. We have seen in Eq. 
(|21|l that these fluctuations result in a phase resolution 



Ae^ 



(24) 



where we have neglected the shot-to-shot number fluctu- 
ations for s impli city. The density noise is here propor- 
tional to (Gi), allowing at best for interferometry at 
the standard quantum limit. 

Figures El and |31 shows the ground state distributions 
in terms of both the probability amplitudes |cm|, and of 
the amplitudes \ce^ \ of the ground state expressed on a 
basis of so-called relative phase states. 
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and the discrete relative phases are given by 
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(27) 



for an arbitrary reference phase chosen here to be zero. 

For weak interactions, i.e. small |G|, the system is in 
a state reminiscent of a coherent state, with a relatively 
well defined phase with fluctuations consistent with the 
standard quantum limit. For large attractive interactions 
the ground state approaches a double peaked distribu- 
tion, corresponding to a macroscopic superposition state 
as described in |2a, |23| . For G ^ — 1 , the ground state 
can thus be approximated by 
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This gives for the phase distribution 
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FIG. 2: Ground-state distributions as function of the dimen- 
sionless parameter G for attractive interaction. In a) mag- 
nitudes of ground state components |cm(G)| are shown for 
attractive interaction, and in b) the distribution of relative 
phase components lce„(G)| in Eg. 1)26^. 
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FIG. 3: Ground-state distributions as function of the dimen- 
sionless parameter G for repulsive interaction. In a) mag- 
nitudes of ground state components |cm(G)| are shown for 
attractive interaction, and in b) the distribution of relative 
phase components \cb^{G)\ in Eg. (1261 . 



The transition between the superfluid and superposition 
regimes takes place just below G = — 1, where the relative 
phase becomes well defined due to squeezing. 

For large repulsive interactions the ground state goes 
through the Mott-insulator transition in a continuous 
manner. Here the number distribution narrows and the 
phase distribution widens until it becomes essentially flat, 
indicating a completely random phase when averaged 
over an ensemble. 

Figure01shows the uncertainties in the relative number 
Aj|, and the relative phase AJ^, together with their 

product AJyAJz for N = 100 atoms. The uncertainty 
product is essentially constant in the region 

-l<G<iiV2, (30) 

where the upper limit was found from the simulations. 
In this interval the ground state resembles a coherent 
state, making this the superfluid regime. For large at- 
tractive interactions the uncertainty product increases at 
the superfluid-superposition transition and saturates for 
strong interactions. As the uncertainty in relative phase 
does not increase when compared to its value in the su- 
perfluid regime, the increase in the uncertainty product 
can be attributed solely to the increase in AJ^. The fact 



that the ground state in the superposition regime is dou- 
ble peaked, as can be seen in Fig. [21 and thus no longer 
a minimum uncertainty state, gives using Ea. (|28|l for the 
uncertainty in number difference 

{J'.) -J'-^, (31) 

in agreement with the numerical results, see for instance 
Fig. 21 On the Bloch-sphere the distribution is concen- 
trated around the two regions m = ±J, this bimodal 
character giving giving rise to a large spread in AJz- 

Figure|31shows the visibility V of the fringes, Eq. H17|l. 
after the free expansion of the matter waves following 
the switching off of a static double-well potential. In 
case tunnelling dominates the dynamics of the system, 
we observe high-contrast interference fringes, but there is 
an abrupt transition, with contrast decreasing as |G|^^ 
as soon as G < —1. The underlying physics govern- 
ing the decrease of contrast is essentially the same as in 
the superffuid-Mott insulator transition predicted ;23;, 3| 
and observed in optical lattices . If tunnelling dom- 
inates, the state of the two-mode system is essentially 
a superfluid with a well-established phase relationship 
between the two modes of the beam splitter, resulting 
in high-contrast interferences. In the superposition and 
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FIG. 4: The uncertainties AJ^, Ajf and the product 
AJyAJz as functions of G for A'^ — 100. The phase un- 
certainty AJy goes through a global minimum at G = —1, 
the border between attractive superfluid and superposition 
states, and increases in the repulsive regime until it saturates 
to a constant value as the system passes through the Mott- 
insulator transition. 



Mott regimes, by contrast, the two wells are isolated from 
each other, with no phase relationship between them. 

This behavior of the fringe visibility can be further 
understood from the Heisenberg uncertainty relation 



AjyAl^^\{[jy,l 



(32) 



which shows that the uncertainty product AJyAJ^ is 
proportional to the polarization {J^} and thus the vis- 
ibility (|17|l . This shows that the regime of high visibility 
coincides with the one where uncertainty product is con- 
stant. This is, as expected, the superfluid regime. The 
polarization measures the coherence, here the degree to 
which adjacent wave function components Cm are popu- 
lated. In the Mott-insulator and superposition regimes, 
the wave function becomes highly peaked around the cen- 
ter (minimal |m|) and boundary (maximal |m|) which 
thus limits the coherence. 

In addition to illustrating that the visibility decreases 
outside the superfluid regime given by Eq. (|^ . FigEl 
also shows that for G > and even particle number, the 
asymptotic value for the visibility is V = ^ rather than 
zero for N odd. This difference in asymptotic behaviors 
can be explained as follows. For repulsive interactions 
and far into the Mott insulator regime, the ground state 
is approximately be given by a Fock state with equal pop- 
ulations in both wells. In the case of odd atom number, 
the additional atom can be in either of the wells, so that 
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FIG. 5: Visibility of interference fringes for a ballistically ex- 
panding two-mode BEC as a function of G, the ratio between 
the mean-field and tunneling energies, for particle numbers 
ranging from 10 to 200. For attractive interaction the visibil- 
ity decreases abruptly for G < — 1. In the case of repulsive 
interaction the behavior is qualitatively and quantitatively 
different for even (solid) and odd (dashed) atom numbers, as 
explained in the text. 



This gives for the polarization 



(34) 



and hence V = 1/2 as asymptotic value. 

For the case of even particle number, in contrast, the 
asymptotic ground state for G ^ oo is the Fock state 
|m = 0) with N/2 atoms in each well. For large but 
finite values of G the state can be approximated by 



IV) = 71^1 J, 0) 



\J,l) + \J,~l) 



(35) 



where epsilon is a small number. This gives for the po- 
larization 



(36) 



and an asymptotic value of the visibility ^ 0. 

Let us now turn to the phase resolution A8^ of the 
beam splitter. The left-hand side of Fig. shows the 
variance AO^/AGgq^ after free ballistic expansion from 
the two-mode ground state of the beam splitter, plot- 
ted as a function of the ratio G between the mean field 
energy and tunnelling energy. The phase resolution ex- 
hibits a sharp minimum just below G — —1. At this 
point, we find numerically that the energy gap between 
the ground state and the first excited state of the double- 
well goes to zero within our numerical accuracy. This 
indicates the presence of a 'quantum phase transition' — 
or more precisely cross-over — a property also inferred 
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FIG. 6: Squeezing of the phase variance relative to the stan- 
dard quantum limit for the two-mode ground state. For at- 
tractive interaction the phase squeezing is minimal just below 
G = — 1. The crossover from Poissonian fluctuations in the 
superfluid regime G < — 1 to the superposition state G > — 1 
becomes sharper with increased number of particles. For re- 
pulsive interaction there is a smooth transition between the 
superfluid regime of low G- values to the Mott-insulator state. 
In the asymptotic limit G ^ oo, the squeezing approaches a 
plateau, which is different for odd and even atom numbers. 



from Fig. |2] where the quantum statistical properties of 
the ground state obviously change at the point G = — 1. 
This transition is associated with the onset of increased 
number fluctuations between the populations of the two 
modes, and correspondingly to suppressed phase fluctu- 
ations. Hence, the phase resolution of the beam splitter 
is maximized at that point. As expected, the crossover 
from Poissonian fluctuations in the superfluid regime for 
G > —1 to a superposition state for G < — 1 becomes 
sharper with increased particle number. For attractive 
interaction there is thus a squeezed state at the bound- 
ary between the superfluid and superposition regimes. 
For small atom numbers and attractive interaction, the 
minimum of AO^ occurs to the right of the superfluid- 
superposition transition. From this discussion, it might 
appear favorable to operate the beam splitter in that 
regime. This is however misleading, as we must also take 
into account the fact that the fringe visibility rapidly de- 
creases in that regime, as shown in Fig. [S] The situation 
improves rapidly for large iV, though. The sharpness of 
the minimum, however, makes this state challenging to 
create experimentally in a controlled manner. In addi- 
tion, the maximum achievable squeezing is limited due 
to the size limits imposed by the metastability of attrac- 
tive condensates. 

The situation is slightly more subtle in the case of re- 
pulsive interactions. From Fig. 2, it is quite clear that 
as the system moves into the Mott regime, the phase un- 



FIG. 7: Phase variance in the ground state of the double- well 
system at the "critical point" G = — 1 for attractive conden- 
sates, and at G = for repulsive condensates, as functions 
of particle number. The dashed lines show the standard quan- 
tum limit (SQL) and Heisenberg limit (HL) for comparison. 



certainty increases and the number fluctuations of the 
ground state become more strongly squeezed. This is 
also evidenced in Fig. 3, which shows the monotonic in- 
crease in AJy for increasing G. Hence, it would appear 
that the phase resolution becomes increasingly worse in 
this regime. This difficulty can however be eliminated by 
using quantum tunnelling to turn number squeezing into 
phase squeezing, as discussed in section II. As a result, 
the minimal obtainable phase fluctuations for G > are 
given by 



AJ| 



(37) 



since AJ^ — (J^) is the minimal uncertainty of in 
deep the Mott regime. The right-hand side of Fig. El 
shows that the phase resolution achieved by this tech- 
nique in the case of repulsive interactions. Just as for 
the visibility, we find here different asymptotic behavior 
for odd and even particle numbers with a factor of two 
in difference for large interactions. Using Ea. H33fl we find 
for the uncertainty in number difference 



iJi) = 
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which together with Ea. (l34|l for the polarization gives 
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in agreement with the asymptotic values for odd particle 
numbers, shown dashed in Fig. |^ For even numbers, 



8 



0.75 



0.5 



0.25 



T = 500 




1000 



2000 3000 4000 5000 

Time 



FIG. 8: Phase squeezing for g < when AE{t) is taken 
to vary according to Eg. 14211 . for several values of the r. The 
dashed lines show squeezing of the instantaneous ground state 
for comparison. 



Eq.(|2nil yields 
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which together with Ea. (|36|l gives 
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in agreement with Fig|Sl The behavior of both phase res- 
olution and visibility for odd versus even particle num- 
bers can thus be understood from the form of the ground 
state in the limit of large repulsive interactions. 

Strictly speaking, AJj, is the phase uncertainty only 
when the distribution spans an area narrower than the di- 
ameter of the Bloch sphere. For larger values the spread 
A Jy behaves differently, as it is the projection onto the 
plane spanned by Jy and Jj, and not the distance along 
the equator. For the interferometric setup considered 
here the phase uncertainty is in general not equal to the 
phase resolution given by Eq. H21I) , as the latter will grow 
when the visibility goes down. 

We have seen in Section II that the phase sensitivity 
of the system is closely related to spin squeezing. The 
squeezing of both optical and matter waves is of consid- 
erable interest in interferometry, as it offers the potential 
to beat the standard quantum limit of detection, and pos- 
sibly replace it by the so-called Heisenberg limit scaling 
as Fig. [3 shows the phase variation of the ground 

state at both the "critical point" G = — 1 for attractive 
interaction, and at G = \N'^ for repulsive interaction 
as functions of particle number N , the dashed lines giv- 
ing the standard quantum limit (SQL) scaling as l/ViV 
and Heisenberg limit (HL) for comparison. The G- value 
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FIG. 9: Squeezing of the relative phase AO for a double well 
split in time according to equation Eq. 1421 . 



in the repulsive regime was chosen to produce maximal 
squeezing for odd particle numbers, slightly worse than 
for even numbers. Since the particle number can not be 
controlled down to single units, this produces a conser- 
vative limit for the achievable squeezing. 

The combination of the mean-field interaction and tun- 
neling clearly leads to the squeezing of the matter-wave 
field and a sensitivity significantly improved from the 
standard quantum limit under the conditions mentioned 
above. We find numerically that for attractive interaction 
A8^ scales as N~^-^^, slightly above the asymptotic limit 
of N~^/^ derived by Kitagawa and Ueda |[31j. The differ- 
ence is expected to decrease with increased particle num- 
ber. For repulsive interaction we find A8^ oc A^~^, thus 
corresponding to Heisenberg limited phase resolution. 
This comes at the price of larger interactions and also 
necessitates an additional rotation of the Bloch sphere. 



B. Dynamics 

The previous section investigated the ground state 
properties of the static double-well system. Here we 
discuss ways to create states with the desired proper- 
ties starting from a condensate in the ground state with 
a small interwell separation. The adiabatic theorem of 
quantum mechanics states that a system governed by 
a time dependent Hamiltonian and initially prepared in 
an eigenstate will remain in the instantaneous eigenstate 
given that the Hamiltonian changes sufficiently slowly. 
The aim here is evolve the system adiabatically '33^ into 
states of maximal squeezing and then freeze the dynam- 
ics. 

Controlling the magnitude and sign of two-body inter- 
action is readily achievable using Feshbach resonances. 
The rate at which the magnetic field can be swept across 
a resonance is much higher than typical trap frequencies 
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and tunnelling rates, allowing for a practically instanta- 
neous switching of the nonlinear interaction. Decreasing 
the tunnelling rate can be achieved by rapidly separating 
the two wells as the point of minimal phase fluctuations 
is reached. 

We model the splitting of the symmetric double-well 
potential using a tunnelling potential AE{t) that de- 
creases exponentially in time 

AE{t) = hu exp - ^ - ri j , (42) 

where dmm is the minimal separation between the wells, 
and the constant F is a measure of how fast the wells 
are separated. The separation is assumed to occur over a 
finite time r and to reaching a value d(r) such that G ~ 
— 1 for attractive interaction, and G = ^N"^ for repulsive 
interaction. The constant F is adjusted accordingly. As 
the time t = r is reached, the two-body interaction is 
rapidly taken to zero to freeze the dynamics. 

The evolution after t = t is linear and solely governed 
by tunnelling. Note that for attractive interactions, the 
residual rotation of the distribution resulting from quan- 
tum tunnelling has to be kept small in order to avoid 
transforming the phase squeezing into number squeezing, 

/C30 
AE{t)dt < hi:, (43) 

This can be achieved by either taking A.E{t) ^ 1 at the 
point where G = — 1, or separating the wells rapidly after 
i = r to make the left-hand side of Eq. H43|l negligible. 

The squeezing in the phase variation for attractive in- 
teraction is shown versus time in Fig. |Slfor several val- 
ues of the parameter t. The instantaneous squeezing of 
the true ground state is shown dashed for comparison. 
The evolution of the system is initially adiabatic, but 
becomes diabatic for larger separations making the sys- 
tem freeze out in a state less squeezed than the desired 
ground state. Longer evolution times make the system 
more adiabatic thus achieving squeezing closer to, but 
still far from the optimal values. For all cases shown in 
Fig. |S| it is however apparent that the dynamics is only 
partially adiabatic and that the minimal phase variation 
is severely limited by violation of the requirement that 
the dynamics be adiabatic. 

Figure|51shows the squeezing dynamics of a condensate 
with repulsive interaction when the tunnelling energy is 
decreased exponentially in time according to Ea. (|42|l . As 
the value G = ^N^ is reached, the two-body interaction 
is rapidly tuned to zero. After that a pulse of tunnelling 
interaction, which corresponds to bringing the two wells 
closer again, is applied in such a way that 

AE{t)dt^—, (44) 

thereby rotating the Bloch-sphere distribution by an an- 
gle 7r/2 around the J^j-axis and thus transforming the 



number squeezing into phase squeezing without other- 
wise changing the distributions. 

Initially the phase variance is seen to increase as the 
state becomes number-squeezed for increasing values of 
G{t) until the time i = r is reached when g ^ 0, ef- 
fectively freezing the evolution until the tunnelling pulse 
(|44|l is applied. For all values of t, which are indicated in 
FigEl the same tunnelling pulse was used. The dynamics 
is seen in Fig. |51to become more and more adiabatic as 
the value of r is increased, just as in the case of attrac- 
tive interaction. Here, however, the maximal squeezing 
achieved is larger than for the attractive case. This is not 
obvious at first, even though the repulsive ground state is 
more squeezed, as the system has to be evolved to values 
of stronger interaction than in the attractive case. 

IV. CONCLUSION AND OUTLOOK 

In this paper we have investigated the full quantum dy- 
namics of a condensate in a symmetrically split double 
well. Expressions for the visibility and phase resolution 
during a ballistic expansion stage were given and inves- 
tigated numerically. The possibility of creating phase 
squeezed ground states by adiabatic splitting was demon- 
strated, but is limited by the time scales involved. The 
increased phase sensitivity for ground states of attractive 
condensates, which are known to be stable only for parti- 
cle numbers [2^ up to around A'' « 10"^, was found to re- 
quire long splitting times to achieve adiabatic evolution. 
For repulsive condensates where the densities are limited 
only by the requirement that the two-mode model be ap- 
plicable, a scheme with scaling at the Hcisenberg limit 
was outlined and tested in dynamics simulations. The 
original Heisenberg-limited scheme suggested in 16] for 
optical MZ interferometers and later applied in various 
forms to the case of atomic condensates llE ' ^s^'^ 
dual Fock-states. i.e. the state was of the form given 
by Ea. l|5S|l . In the present context such states are un- 
suitable due to their low visibility. The problems can be 
avoided by not taking the system deep into the Mott- 
insulator regime where the ground state is the dual Fock 
state, but rather to the threshold where G — N'^/2, and 
the visibility still is high. 

It is well known that the introduction of a linear poten- 
tial, for instance due to gravity, changes the localization 
properties of the double-well eigenfunctions. The corre- 
sponding effects in the many-body regime are presently 
explored experimentally 0| and will also be the subject 
of future theoretical investigations. 
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